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Abstract 

Using the SU (5) symmetry of the J^D hyperdiamond and results on the study 
of 4D graphene given in |10| . we engineer a class of J^D lattice QCD fermions whose 
Dirac operators have two zero modes. We show that generally the zero modes of the 
Dirac operator in hyperdiamond fermions are captured by a tensor fl l with 4x5 
complex components linking the Euclidean SO(4) vector ju; and the 5-dimensional 
representation of SU (5). The Borigi-Creutz (BC) and the Karsten-Wilzeck (KW) 
models as well as their Dirac zero modes are rederived as particular realizations of 
0^ . Other features are also given. 

Keywords: Lattice QCD, Borigi-Creutz and Karsten-Wilzeck models, J^D hyper- 
diamond, Jj.D graphene, SU (5) Symmetry. 



1 Introduction 

One of the highly important things in lattice QCD is the need to have a fermion action 
with a Dirac operator T> having two zero modes at the points K and K' of the reciprocal 
space; so that they could be interpreted as the two light quarks, up and down. In this re- 
gards, Michael Creutz proposed few years ago an interesting J^D lattice fermion action [T] 
by extending the dispersion energy relation of 2D graphene [2j [3] to J^D quaternions. The 
Creutz fermions, which exhibit very useful properties for dealing with lattice QCD, was 
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further developed by Artan Borigi [I] leading afterwards to the so called Borigi-Creutz 
(BC) fermions [6]. This is a simple J^D fermion lattice model using a 4-component 
Dirac spinor \l/ r as well as a particular linear combination of the 4x4 matrices 7„. BC 
lattice model has too particularly chiral symmetry, preventing mass renormalization, and 
a Dirac operator V BC with two zero modes K BC and K' BC that describe a degenerate 
doublet of quarks. The engineering of these two zeros is mainly due to mimicking 2D 
graphene as well as to a nice trick relying on using a " complexification" of the Dirac 
matrices type T M = 7^ + i^'^ and the remarkable linear combination 7i + 7 2 + 7 3 + 7 4 - 
This way of doing let suspect that the BC model hides a more fundamental property 
that controls, amongst others, the engineering of Dirac operator zeros. 

In this paper, we contribute to this matter by proposing a fermionic 4D lattice QCD-like 
model based on SU (5) symmetry and containing BC and KW [7J |8] fermions as partic- 
ular cases. In our proposal the J±D lattice, denoted as L 4 , is given by the hyperdiamond 
generated by the 4 simple roots a%, a 2 , a 3 , a 4 of SU (5) [9]. Moreover, each site of L 4 
has 5 first nearest neighbors (r^ + aXi) forming a 5-dimensional representation of the 
SU (5) symmetry. 

The proposed free massless fermion lattice action is given by 

4/4 4 \ 

s - 4 EE E*' W) +E*^ (f «i) ( L1 ) 

r ,u=l \l=0 1=0 J 

where \l/ r lives on L 4 and where the 5 vectors A , Ai, A 2 , A 3 , A 4 stand for the weight 
vectors of the 5- dimensional representation of SU (5) satisfying the basic property A + 
Ai + A 2 + A3 + A 4 = 0. This identity has a physical interpretation in terms of conservation 
of total momenta at each site r; and has as well a strong link with the linear combination 
7i+72+73+74~ 2r = used in BC model. The distance a ||Aj|| = a^0- is the parameter 
d of the lattice. 

Notice that our lattice fermion model depends on the complex tensor Q l whose features 
will be studied here in details. It is a kind of a complex 4x5 rectangular matrix that 
plays a quite similar role to the matrix combinations 7^ + ij' , 7^ + «7 4 (1 — # m4 ) and 
7i + 72 + 73 + 74 used in the BC and KW models [5J E]. In the case of free BC fermions, 
we will show that the corresponding (f2jj BC tensor reads as 



/ 2 l + 2i -1 -1 -1 \ 

2 -1 l + 2z -1 -1 

2 -1 -1 l + 2i -1 

\ 2 -1 -1 -1 l+2i J 



;i-2) 



2 



and for the KW ones, the associated (fijjj RW is given by 

/ i \ 

i 

i 

\ 2 1 1 -2 z - 2 / 



1.3) 



Notice also that, generally speaking, Q 1 ^ transforms in the bi-fundamental of SO (4) x 
SU (5); and should not be a tensor product; i.e ft! ^ £^ x ( l in order to recover the 
continuum limit. This Q 1 ^ links then the 5*0 (4) and SU (5) fundamental representations 
respectively captured by the indices /i = 1, 2, 3, 4 and I = 0,1, 2, 3, 4. It allows as well to 
implement quark-quark-gluon interactions by letting it to depend on positions r. 
The tensor Q 1 ^ has another remarkable feature that we are interested in here; it plays 
a crucial role in the engineering of lattice Dirac operator with a generic number N of 
zeros; including the interesting case N = 2. Splitting Q 1 ^ into real and imaginary parts 
as R l + i J 1 , it is clear that the set of real parameters captured by the rectangular matrix 
B; , and similarly for J 1 , is given by the real 20-dimensional moduli space 

m _ SO (4, 5) 

^ 20 " SO (4) x SO (5)' (L4) 

so that the total moduli space of fll.ip is given by the complexification of 9^ 20 namely 
u(i)xu(5) wr th dimension diniR U (4, 5) — dim# U (4) — diniR U (5) = 40. As such, the zero 
modes of the lattice "Dirac" operator V^q ~ 7 M (f^e mk ' Ai +fi^e~ mk A! ) of the above 
lattice action (11.11) correspond to some particular tensors (Q^) ; and so would live 
on subspaces of 9^20 x ^20 • From this moduli space view, operators V^n that have two 
zero modes have tensors type (^) 2zeros which require the fixing (4 + 4) real moduli; 
and therefore would live on real 32-dimensional subspaces 9JT32 of the moduli subspace. 
Typical examples of 9Jt 32 is given by the particular space 

U ^ (1.5) 



U(A) x C/(4) 



or more generally by x with n + m = 8. In addition to the previous BC and 

KW examples with two zeros, one can engineer others by solving appropriate constraint 
relations or directly by making wise choices based using symmetries. This analysis is 
developed in sub-section 4.2 from which we learn that the real (B 1 ^ BC and imaginary 
W) bc P ar ^ s °f the BC complex tensor (f^) should obey the constraint relations 



J2( r »)bc = -(K)bc > E(4 C -(* C 
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Similar relations can be also written down for the tensors (R l „) „„, and (Jj.) _,„, in the 
case of the ifW model; see eqs f)4.45j) . 

The presentation of this paper is as follows: In section 2, we recall briefly some results 
on the BC and KW models; in particular the engineering of the two zeros of the cor- 
responding Dirac operator. We also give extensions of these models where the linear 
combination r 5 = r y l + 7 2 + 7 3 + 7 4 is interpreted in terms of a fifth direction; the four 
others are T M = 7^ + ij' . These extensions may also be understood as an indication for 
the existence of the tensor Q l that captures data on the zeros of Dirac operators. In 
section 3, we develop our proposal; and show, amongst others, that the Dirac operator 
following from (11.11) can also have two zero modes by making wise choices of Q l . In 
section 4, we re-derive the BC and KW theories, and their extensions using Ts term, 
from ( 11. ip . In section 5, we give a conclusion and a comment on gauged ( 11.11) . 



2 BC/KW fermions and beyond 

First we study the two zero modes Kbc and K' BC of the Dirac operator T>bc of the BC 
fermions; then we give those K KW and K' KW of the Dirac operator V KW of KW model. 
We also use this section to give an extension of BC and KW models where the term 
* r r^ r is modified as ^> r T^ r+ab . 

The analysis given in this section is useful to fix the ideas; but also for later use since 
BC, KW models and their extensions can be re-derived from the lattice action (11. ip by 
using wise choices of the tensor £l l . 

2.1 BC fermions: the free model 

Following jH |6] and using the 4-component Dirac spinors = (0">Xr)> the lattice 
action of free BC fermions reads in the position space, by dropping mass term m , as 
follows: 

Sbc - (E - £ f r+aM T^ r J - f £ * r T* r ^ 



where, for simplicity, we have dropped out gauge interactions; and where and T M = 

7 M + 27'^; which is a kind of complexification of the Dirac matrices. 

Moreover, the matrix T appearing in the last term is a 4 x 4 matrix linked to 7 M , 7 //x as 
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follows: 

y*» = T - 7^ , 7^7" + 7^7^ = 26'"' 

4 

2T =^2Y , 7 M + «7 ,M = 
Mapping (12.11) to the reciprocal space, we have 



(2.2) 



Sbc ~ ^ ^k^Bc^k (2.3) 

k 

where the massless Dirac operator T>bc is given by 

+^(^ + Tjsin(a^)-fr . 

Upon using T /t + = 27^ and Y^ — = 2ry^, we can be put Vbc in the form 

V BC = D^ + D u --T (2.5) 
a 



(2.4) 



with 



D * = « ( $^7 M sinafcJ , D k = if £ cos afc„ 

vM =1 / \M =1 



4 

V^cosrxfc,. I 

(2.6) 



where A; M = k./x. In the next subsection and in section 5, we will derive the explicit 
expression of these k^s in terms of the weight vectors A; of the 5-dimensional represen- 
tation of the SU (5) symmetry as well as useful relations. 

The zero modes of Vbc are points in the reciprocal space; they are obtained by solving 
Vbc = 0; which leads to the following condition 

Y (sin aK^ - cos aK^) - V [2 - ^ cos aK n ] = ■ ( 2 - 7 ) 

This condition is a constraint relation on the wave vector components K^, it is solved 
by the two following wave vectors: 

point K BC : K 1 = K 2 = K 3 = K 4 = , 
pointy : K[ = K> 2 = K>z = K>± = f a , 

and are interpreted in lattice QCD as associated with the light quarks up and down. 
Notice that if giving up the 7^- terms in eqs fl2.m2T3]) ; i.e 7^ — > 0, the remaining terms 
in Vbc namely ~ 7^ sin aK^ have 16 zero modes given by the wave components 
= 0, 7r. By switching on the 7^-terms, -/^ zeros are removed. 
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2.2 Beyond BC model 

Although very important, the lattice action Sbc is somehow very particular; it let sus- 
pecting to hide a more fundamental property which can be explicitly exhibited by using 
symmetries. This will be developed in details in next sections; but in due time, notice the 
three following features that can be understood as a way to exhibit the SU (5) symmetry 
in BC model. 

1) The price to pay for getting a Dirac operator with two zero modes is the involvement 
of the complexified Dirac matrices T M , as well as the particular matrix V . 
Despite that it breaks the SO (4) Lorentz symmetry since it can be written as 



v, 



(2.9) 



with 



v. 



( 1 \ 

1 

1 

V i / 



(2.10) 



violating explicitly the SO (4), the matrix T plays an important role in studying 
the zero modes. The expression of the matrix F (12. 2p should be thought of as 
associated precisely with the solution of the constraint relation 

4 

2r-^ 7 M = (2.11) 

that is required by a hidden symmetry of the BC model namely the SU(5) sym- 
metry of the 4D hyperdiamond; see also the correspondence given in section 4 of 
ref HDI. 



2) the BC action Sbc lives on a 4D lattice Lf c generated by /x = v ; i.e the vectors 



( vS \ 



( vS \ 



( v x 4 \ 



Vi= 



v 2 = 



v 3 = 



v 4 = 



(2.12) 



\vlJ 



\vlJ 



These /x- vectors look somehow ambiguous to interpret by using the analogy with 
4D graphene prototype; they are also ambiguous to interpret from the SU(5) sym- 
metry view. Indeed, to each site 
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there should be 5 first nearest neighbors that are rotated by SU (5) symmetry. 
But from the BC action we learn that the first nearest neighbors to each site r are: 

r+av x 
r+avo 

r -> < , (2.13) 

r+av 3 

k r+av 4 
The fifth missing one, namely 

r -)> r+av 5 , (2-14) 

may be interpreted in the BC fermions as associated with the extra term involving 
the matrix V. 

3) To take into account the five nearest neighbors, we then have to use the rigorous 
correspondence 

->■ v M , and T 5 -> v 5 (2.15) 
which can be also written in a combined form as follows 

r M ^v M , (2.16) 

with 

r M = (F\r 5 ) and v M = (v At ,v 5 ). (2.17) 
Because of the .St/ (5) symmetry properties [9], we also have to require the condition 

Vi + v 2 + v 3 + v 4 + v 5 = , (2.18) 

characterizing the 5 first nearest neighbors. 

To determine the explicit expressions of the matrices Tm in terms of the usual 
Dirac ones, we modify the BC action (12. 1 p as follows 



5 5 

v w r M y„ \ . 

(2.19) 



r \Af=l M=l 

exhibiting both SO (4) and SU(5) symmetries and leading to the following free 
Dirac operator 

4 

V = Ta H ( r M + r J sin (akfj) + ± (F 5 + f 5 ) sin (ak 5 ) 

(2-20) 

^2 ~ *V) cos (°^) + i ( r 5 - f 5 ) cos (ak 5 ) 



2a 



where ku = k.v^ and where 

5 5 
M=l M=l 

expressing the conservation of total momenta at each lattice site. 
Equating with ( 12.4112. 51123]) . we get the identities 



t^ ^ /j, ~\~ , 
t — t = r — f 

and 

i a (r 5 + f 5 ) sin {ah) + £ (r 5 - f 5 ) cos {ah) = -fT 



T 5 = — 2iT for sin {ah) = 
T 5 = —2T for sin (ak 5 ) = 1 

where k 5 = — {h + k 2 + k 3 + h). 



(2.22) 



(2.23) 



Eqsf l2.22j) are solved by T M = Y M ; that is 

T M = 7 M + *7' M = Y + * (T - 7 M ) (2.24) 

while 



(2.25) 



In this 5-dimensional approach, to be further be developed in next sections, the ambi- 
guity in dealing with the ft- vectors is overcome; and the underlying SO (4) and SU (5) 
symmetries of the model in reciprocal space are explicitly exhibited. 

2.3 KW fermions 

The BC fermions described above is not the only model that has a Dirac operator with 
two zero modes. There is an other remarkable model that has been considered in recent 
literature [5], [6] . This is given by the so called Karsten- Wilzeck model [7J [8] whose free 
fermion lattice action reads as follows, 

S KW ~ ( E ^ r W ~ ^4 (1 - Ml ^ r+a A + £ * r [^0 + f 7 J 

r \/x=l / r 

-S E (E *W [7* + *7 4 (1 - Ml ^ • 

(2.26) 
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This is quite similar to (12. ip ; it is obtained by making the following substitutions 

y M _ ► 74 (i - ^ , r — ► f 74 ■ ( 2 - 27 ) 

Moreover mapping the above lattice action to the reciprocal space, the Dirac operator 
T>kw reads as, upon dropping the bare mass m , as follows: 

4 3 

V KW = sin ak tl + ^7 4 ^(1 -cosa^) = 0. (2.28) 

/x=i fi=i 

This operator has also two zero modes located at the points Kkw and K' KW , 

K KW =(0,0,0,0) 
Kkw = (0,0,0,*) 

The comments we gave in the case of BC fermions apply as well to this case. In particular, 
we have the extended-KW action 

s 'kw~ ^^(^^ r r M vi/ r+aVM -^^ r+aVM r^J , (2.30) 

r \M=1 M=l J 

where now 

r„ = 7 M + ^74 (i - M 

r 5 = -f 7 4 > for sin (a/c 5 ) = , 
F 5 = ~f74> for sin (aA; 5 ) = 1 . 



In section 5, we will show that the BC and KW fermions as well as their extended 
versions with lattice actions S' BC and S' KW respectively given by (I2.19p - fl2.30p . correspond 
in fact to particular realizations of the tensor Q l of the proposal (II. ip . 



3 Hyperdiamond fermions 

Before describing the proposal ( 11. ip . it is interesting to start by recalling useful features 
of the hyperdiamond L4 [9j QUI E] . 



3.1 the 4D lattice 

The lattice L 4 is a J^D lattice made by the superposition of two sublattices A and B with 
respective coordinates positions, 

r -f= Ntct! + N 2 cx 2 + N 3 a 3 + N±a 4 , 
r B = r A + d^f\ t , a = d& 

9 



where the iVj's are arbitrary integers and the A/'s are the weight vectors of the 5- 
dimensional representation of SU (5). For a related construction see [T5l ITS] . 
Recall that each A; weight has 4 components 

Xf=(Xl,XlXf,Xt), 1 = 0,1, 2, 3, 4, (3.2) 

obeying the group representation identities 

A£ + A? + A£ + A£ + A£ = , / u = l,...,4 . (3.3) 

These identities can be interpreted physically in terms of the conservation of total mo- 
menta at each site of L4. 

Ilexp (f X>X) = fl ex P (f X>X) = 1 • ( 3 - 4 ) 

Notice that the components Af can be expressed in various, but equivalent, ways. Ex- 
plicit ways were used in [91 [TOj [TTJ; see also [12], but, to our understanding, the more 
convenient and significant way to do is the one using the SU (5) simple roots as follows: 

A = + §Q!2 + |a 3 + |a 4 



Ai = -iai + |a 2 1 L '"' 1 
+ §«3 + |c 

A 4 = - §t*2 - §0:3 - f £*4 



5 1 1 5 z 1 5 J 1 5 
|ai - f a 2 + §a 3 + §c~. 



A 2 = - 5«2 + 5CK3 + 5«4 

A 3 = -iai - fa 2 - |a 3 + 



These vectors capture the fact that each site in the sublattice A has 5 first nearest 
neighbors in the sublattice B, 

one site — > 5 first nearest neighbors 







r A + aX 






r A + aX 1 


r A eA 


< 


r A + aX 2 G i3 






r A + aA 3 






r A + aA 4 



(3.6) 



Notice also that the expression of the A;'s in terms of the roots is helpful not only for 
calculations; but also for exhibiting manifestly the role of the SU (5) symmetry. 
For later use, we also need the dual vectors, 



^ 27^ 2nN 2 2ttJV 3 2nN A 
G = u>i H w 2 H ^3 H W4, 



(3.7) 
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obeying the reciprocal lattice property 

exp {iaG.r A ) = 1, 
thanks to the duality relation uJi.ctj = 5ij. 



(3. 



3.2 the free model 



The proposal (11.11) involves, in addition to the Dirac spinor \I/ r and the 7^ matrices, the 
remarkable complex tensor Q l • this is a bi-fundamental representation of SO (4) x SU (5) 
given by the complex 4x5 matrix, 



/ n° n\ nj nf nf \ 

$~^2 ^"^2 ^~^2 ^~^2 ^^2 

^"^2 ^"^3 ^~^3 ^"^3 

V n\ Ql Ql nj J 



(3.9) 



that allows to link the lattice Euclidean space time index /x, capturing a 4-component 
vector of 5*0(4); and the index I of the 5-dimensional representation of the SU (5) 
symmetry of the hyperdiamond. In some sense, it plays the role of 7^ and T in the BC 
model. 

Notice that the action (II. ip can be also written in a similar form ( I2.19p . as follows 

/ 4 4 

S ~ fa E E ^ r r^ r+aA; + £ 9r+»F% 



(3.10) 



where now the T l has the following general form 



E^ 



(3.11) 



v^t=l / \A 1 =1 

In the reciprocal space, our massless free fermion lattice action reads as 

ElE*^ 



s 



(3.12) 



vM=l 



with Dirac operator given by 



(3.13) 



and where and its complex adjoint read as follows: 

4 4 

11 



iak.A 



(3.14) 



So the zero modes of the Dirac operator (13.131) are obtained by solving the following 
constraint equation 



Re XX e H =0 > 



.1=0 



which reads explicitly like 



E K + fy] cos ^ +*K- fy] sin ^ = °> 



(3.15) 



(3.16) 



1=0 



where are the relative phases of the wave functions ^ T +aX t with respect to \l/ r ; they 
read in the first Brillouin zone as follows 



a(k.A,) , 0<$ ; <27r, 



with 



(3.17) 
(3.18) 



k.A, = X> m .a? ■ 

11=1 

to be denoted as fcj = k.Aj. Moreover, because of the SU (5) symmetry (13.31) . we also 
have the constraint relation 

4 

^$, = 0, mod27r (3.19) 

i=0 

showing that the 5 phases are related; so that one of them; say $ , can be expressed 
in terms of the others like 



$o = - ($i + $2 + $ 3 + $4) • 
Notice that if solving the constraint eqs ( 13.161) as 

$, = ^, = 0,1,2,3,4, 



(3.20) 



(3.21) 



(3.22) 



we end with a constraint relation on Q l namely 

I>(ny = o. 

1=0 

This condition can be solved by using eq (13.3j) ; this shows that Re (^1) can be thought 
of as proportional to the weight vectors Xf of SU(5). A particular realization of is 
given by the one considered in [9J [TUl [11] ; and which reads as follows 



( 





V5 


v/5 


-V5 


-V5 \ 







v/5 


-V5 


-V5 









v/5 


-v/5 


v/5 


-V5 


\ 


-4 


1 


1 


1 


1 / 



(3.23) 
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Notice in passing that the tensor (fiJJ and (f^J RW given in the introduction, and 
which are associated with the BC and KW models, obey as well the property (13.22)) : 
this feature will be proved in next section. 

Below, we want to look for general solutions going beyond $; = ^^; and recover the 
special solutions ( I3.2ip as particular cases. This also allows us to get more information 
on the role of the tensor f2„. 

(J, 

4 Solutions with two zero modes 

First we study a simple model with two zero modes having similar properties as the BC 
and KW ones; then we analyze the generic case involving the 40 real moduli captured 
by the complex tensor f2^; and explore the link between the number of zero modes and 
the moduli space of Q 1 ^. 

4.1 Simple model 

This lattice model is a simple toy prototype having a lattice Dirac operator T>(k), of 
quite similar form to BC and KW operators considered in section 2, sharing with eq fll.ll) 
the three following features: 

• the operator T> (k) has two zero modes located at points K = (Ki, K2, K%, K4) and 
K'= (K[, K' 2 , K' 3 , Ki) of the reciprocal space; that is: V (K) = V (K 7 ) = 0, 

• the continuum limit of T> (k) in the neighborhood of the points K and K' is given 
by the usual free Dirac operator ^ M 7 M ^ plus higher order terms that have an 
interpretation in the Symanzik effective theory. 

• it has the SU (5) symmetry of the hyperdiamond that permits to mimick tight 
binding model of 2D graphene; there the underlying 2D lattice is given by the 
honeycomb which is known to have an SU (3) symmetry. The SU (5) encountered 
in the present study is just the extension of the SU(3) to J^D dimensions. 



The lattice action S of this model is of the form (11. lft ; but with a very particular tensor 
Q . It reads in the reciprocal space like 

S= j 7^4* (k) V (k) * (k) . (4.1) 
J (2tt) 
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For simplicity, we often write this action like ^k^k^k with the matrix operator "Z\ 
given by the following 4x4 matrix, 



sin ak„ 



— [ cos ak u — v a cos( ak u ) 



(4.2) 



with the vector = (1, 1, 1, 1) as in eq (12.10|) . In addition to the four gamma matrices 
7 M , this operator involves also the particular combination ^^J^v^ required by SU (5) 
symmetry of the hyperdiamond. By comparing this matrix operator with the generic 
one (11 .ip depending on the complex tensor Q l = R l + iJ 1 ^ namely, 

4 



V 



sin aki — R 1 ^ cos aki 



1=0 



1=0 



(4.3) 



we find that (14. 2p is indeed a particular operator of the general (I4.3p . The corresponding 
matrices are as follows, 



R 



m 1+V2 



( 1 -1 \ 
10-10 

10 0-10 

V 1 -1 / 



/ 1 \ 

10 

1 

V 1 / 



(4.4) 



Notice that these matrices have no free degrees of freedom, since all parameters are fixed 
to defined numbers; they correspond to a particular point in the moduli space fj^^^m 
and obey the following remarkable properties 

4 4 

XX = o , £j); = (i + V2)r; • (4.5) 

1=0 v=\ 

These two relations have a physical interpretation; they will be recovered in next sub- 
subsection when we study the zero modes of the matrix operator for generic Q l . 
Eqs (l4.5p are the conditions that should be obeyed by Q l — R l + iJ 1 ^ to have two zero 
modes located at 

(*m) = , «o = ^ 



with 



fj-=l M = l 
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Notice also that for = 0, mod the operator T>k vanishes exactly due to 



(1 + V2) 



tv 



(4.8) 



and for = ^, mod it also vanishes because of the identity 



V2 1 
~2~ ~ 1 + V2 I 2 



(4.9) 



From these relations, one learns that the operator 2\ has, up to translations in the 
reciprocal space of the hyperdiamond, two zero modes located at the points K and K'. 
To see that these wave vectors belong indeed to the first Brillouin zone in the reciprocal 
space 7Z 4 , we first use the expansions 



K = QiVi + Q 2 uj 2 + + 



(4.10) 



with the Qi's and Q"s real numbers; and where the u>j's are the four fundamental weight 
vectors of SU(5) that generate 1Z A ; they obey amogst others eqs (13.7ti3T8|) . Then, we 
compute the K^s and K'^s by help of the the relations ( 13.181) expressing these components 
in terms of the A;'s respectively as K.A^ and K'.A^. Using eqs (l3.5j) giving the A;'s in 
terms of the simple roots a^, ot 2 , 0L3, CK4; and the duality property U)i.ctj 



6ij, we can 



express the components k\ 



( h\ 

k 2 



k.Xi as J2iQi (ui.Xi 
( 



or equivalently in matrix form like, 



(4.11) 



1 

5 


+! 


+1 +i\ 


/ Qi 


\ 


1 

5 


2 
5 


+1 +1 


Q 2 




1 

5 


2 
5 


3 + i 
5 ^5 


Qz 




1 

5 


2 

5 


3 4 / 

5 5 / 


\ Qa 


J 



This leads to 



Qi 
Q 2 
Qs 
Q 4 



(4.12) 



—2k\ — k 2 — k% — &4 
ki - k 2 

k 2 ~ h 
k 3 - k 4 

From these relations we learn the two following: (i) the sum ki + k 2 + k 3 + /c 4 = — k 
is precisely — | (4Qi + 3Q 2 + 2Q 3 + Q4) that follows from the direct use of ( 13. 5ft . (ii) In 
the first Brillouin zone given by — < k^ < -, the two Dirac point K and K' are located 
at 



K = 



K' 



3tt 
4a 



(4.13) 
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Notice that in the neighborhood of K and K', the operator X? q +K and I> q+ K' operators 
behave, up to the first order in the spacing parameter a, as follows 



4 



(4.14) 



2(1+^2) 



where 

v.q = E = Qi + Q2 + Q3 + 94- (4.15) 

v 

These expansions in powers of the spacing parameter of the hyperdiamond give, at zero 
order 0(a°), the usual Dirac operator of the continuum limit. At first order 0(a), 
we have 5-dimensional operators that play an important in the study of the Symanzik 
theory. The latter has a generic effective action that reads in powers of the parameter a 
as follows, 

^//^E^E^ ^ (4-16) 

(Jj 

n j 

Here, the d j) 's are dimension n operators with relevant and marginal ones; 
some of these operators may emerge in the renormalization of the quantum theory. The 
last terms of ( 14. 14ft are linear in a and so lead to O5 operators; for explicit use of these 
operators; see [T2] . 



4.2 Generic solutions 

In this subsection we want to further explore the role of the complex tensor Q l in the 
engineering of zero modes of the matrix operator X?k,f2 of eq fll.il) . We start from eq (!4.3j) : 
then look for the interesting class of those £>k,n operators having: 

• two zero modes located at two different points K and K' of the reciprocal space. 
These modes, which should be non degenerate, require the two set of conditions: 

P K ,n = and D K ',n = , (4.17) 

• the appropriate continuum limit in neighborhood of K and K' given by the Dirac 
operator. From eq ()4.3p . this property requires amongst others that the imaginary 
part of Vt 1 ^ should be as J 1 ^^ A^® B\ 

• the SU (5) symmetry of the hyperdiamond to take advantage of known results on 
2D graphene. 
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To avoid lenghty technical details on the engineering of the zero modes, we will come 
directly to the key points of our method of construction. The basic ideas are illustrated 
through explicit examples given in next sub-subsection; and general results will be given 
later. 



4.2.1 Method of engineering zeros 

First, notice that the values of the wave vector variable k solving = depend on 
the components of Q l ; this feature means that this zero mode equations can be also 
interpreted as constraint eqs relating wave vectors in the reciprocal lattice 1Z 4 to 
points in the real 40- dimensional moduli space 97t 4 o = !tH 2 o x 9^20 given by, 

m gQ(4 ' 5) x gQ(4 ' 5) (AIR) 

Mm - SO (4) x SO (5) SO (4) x SO (5) " ^ 

The relations between the wave vectors k and the moduli ft 1 ., mean as well that we have 
the two following equivalent statements: 

• Given some tensor , corresponding to a point in the moduli space 9JT40, the 
solving of £>k,Q = lead to wave vectors ...j in the reciprocal space 
TZ A . This approach has been used in the above subsection with Q l given by eqs f)4.4p . 

• Inversely, given a wave vector K^, corresponding to point in TZ 4 , the matrix equa- 
tion Pk,si = allow to determine the set of points in VJH40. In other words, Dk,!! 
should be thought of as constraint eqs on the 40 real moduli of the complex tensor 

Below, we develop this second approach. 

Zero mode's constraint eqs and their linearization 
To do so, notice that I^k,q = is a 4 x 4 matrix operator equation that describe a priori 
a system of 16 eqs; but because of the properties of the Dirac matrices, only the 4 of 
these relations are relevant; they are given by, 

Tr (7 M 2> K ,n) = 0, li = 1,2,3,4. (4.19) 

Notice also that in the interesting case of two zero modes of £>k.Q, we need to fix two 
points in the reciprocal space; say K = (Ki, K 2 , K 3 , K±) and K' = (K[, K' 2 , K' 3 , K' 4 ) 
satisying the conditions, 

Tr (7^,0) = 0, Tr (7^,0) = 0. (4.20) 
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These conditions constitute 8 constraints on the free parameters of the complex tensor 
Q l ; their solutions require fixing 8 parameters of Q l in terms of (Ki, K 2 , K s , K 4 ) and 
(K[, K' 2 , K' 3 , K' 4 ); and therefore lead to a reduction of the moduli space WI40 down to 
some subspaces 

9K 4 o -> 9^32- (4.21) 

The mathematical structure of depends obviously on the nature of the constraints 
on Q l . To have an idea on the space 97132, we study below some explicit examples 
where 9JI32 is, for instance, given by the complex coset group manifold ufyx'uii) 01 
the following real one so S (4)xso(3) x so S (4)xSO(5) • ^ eca ^ by the way that dim U (N, M) = 
(N + M) 2 and dim SO (N, M) = \ (N + M) (M + N - 1). 

Substituting the decomposition ft 1 = R l + iJ 1 ^ back into (14.21) . we get the following 
constraint relation for zero modes, 



. 4/4 4 \ 

v ^ = 7^^ E J i sin aki -XX cos ak i = °> ( 4 - 22 ) 

H=l \l=0 1=0 / 

where the aki = &i are the phases of the wave propagation, along the A; direction in the 
reciprocal space, satisfying the SU(5) symmetry condition, 

4 4 

^jfe, = 0, mod *e , ^$ z = 0, mod27r . (4.23) 

1=0 1=0 

Using the diagonal vector v = (1, 1, 1, 1), these constraints allows to express k in terms 
of the others; and similarly for $ - We have: 

4 4 

k = -Y,K = - k - v » $o = -X! $ m = • (4.24) 

fj.=i ^=1 

Notice that eqs f)4.22p is a system 4 equations with 4 variables $1, $2; $35 $4 and real 
moduli given by the composantes of R„ and J'. Moreover because of their dependence on 
the functions cos ^ and sin and also because of the condition f )4.23p . these numerical 
eqs are highly non linear and therefore difficult to solve. To reduce this non linearity, it 
is interesting to work with the following new variables, 

C M = cos^ , S^sin^ , ^ 4 25 ^ 

^ = R^Cq , ip^ = j®s , 

this allow to gain some simplicity into the zero mode constraint eqs. Focusing on the 
first Brillouin zone in the reciprocal space ^ < < ^; and using the above change of 
variables we can put (I4.22p into the form 

4 

^ (i^C, - = JJ5 - i?C , (4.26) 
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but still with the non linear terms 



C-0 — O1O2O3O4 — O1O2O3O4 — O1O2O3O4 — O1O2O3O4 
—01020304 — 01020304 — 0102O3O4 -+- 01020304 

and similarly for others coming from the expansion of So- We will see below how these 
non linearities can be overcome. 
Solving egs ^-lty 

The solutions we are interested in here are those satisfying the 3 properties given in the 
begining of this subsection namely: (i) two zeros modes at K and K' of the reciprocal 
space, (ii) a Dirac operator in the continuum limit of the neighborhood of K and K', 
and (iii) SU(5) symmetry. 

Our method to solve eqs f)4.17l) involves 4 steps as follows: 

(1) Start from the generic zero mode constraint equations T>^ q = given by (I4.26p . 

(2) Take two arbitrary wave vectors K and K' in the first Brillouin zone of the reciprocal 
space, which by help of eqs (l4.25j) . can be related to Ci and Si like, 

e iaK - Xl =Q + iSi , e iaK '- Xl = C[ + iS[ , (4.28) 

leading in turn to the following system of 4+4 relations defining the two zeros 



v f (4-29) 



u=l 



with variables C v and C' v (since S v = ±y/l -C 2 , S' v = ±y/l - C' 2 ) and the free moduli 
R», Rl as well as J£, J°. 

(3) To surround the non linearities, we make a particular choice of and K' in the 
reciprocal space and interpret (14.291) as constraint equations on the moduli. For example, 
one may take the wave vectors and K' as follows, 

^=(0,0,0,0) , ^ . (4.30) 

These vectors are precisely the one that have been obtained in the previous subsection; 
but one can also consider other choices; some of them are described as remarks which 
will be given at the end of this construction. Using the above particular choice, we have 
the following: 

K = , (C ,Sb) = (+l,0) , (C M ,^) = (1,0) 

K'o = ~a ' (^0> ^0) = ( — 1> 0) ' (^V S'n) = O^i ^f) 



19 



(4) Put these values back into eqs fl4.29l) . we obtain the following constraint relations on 
the tensors R\, and Jl, 

4 4 

*?+E^ = °> T, J z = -( 1+ ^ R l ( 4 - 32 ) 

which should be compared with ( 14. 5 p . The first set of constraint eqs reduce the number 
of degrees of freedom of R 1 ^ from 20 down to 16. The second set of constraints do the 
same thing; but with the tensor J^. So the moduli space WI40 = ^20 x ^20 gets reduced 
down to the subspace 

50(4,4) SO (4,4) 



an. 



(4.33) 



32 SO (4) x SO (4) " SO (4) x SO (4) ' 

Therefore, the two zeros (I4.30p live on those missing 8 moduli living on the cosed 
manifold TI40/VJH32 has been freezed by eqsf )4.30l) . 

With this method, one can engineer various 4D lattice models living on hyperdiamond. 
In this regards, let us consider rapidly other examples; and turn after to give a general 
result through three claims. 

A remark and other examples 
The remark concerns the general form of the operator T^n that has two zero modes 
located at the points = (0, 0, 0, 0) and K'^ = (£, £). It is given by: 



irJ2^ [J2 J Z sinak »-J2 R » cosak > 

( 4 ^ 

— J° sin ak.i;+ | ^^-R^ J cosak.t; 



4 



f5> 



K v=l 



with the condition 



E^=0+^)E« 



(4.34) 



(4.35) 



Notice that in the neighborhood of K^, the expansion of £>K+q,n reads, up to second 
order in the power of the lattice spacing parameter a, like: 



V, 



7 f 



t/=l 
' 4 



X)( g ,) 2 i^+(q.t;) 2 i? 



I/=l 



(4.36) 



0(a 2 ) 



where we have used Ylt=i = — • To interpret as the real wave vector of the 
Dirac theory, we have to require moreover 



1 V n, V 70 K U 



(4.37) 
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The other examples concern different kinds of two zero modes other than the ones given 
by (14.301) . For instance, in the case where the two zero modes are chosen like in BC 
model ( 12. 8ft namely 

$1 = $ 2 = $3 = $4 = , $o = , 
^ = $£ = $'3 = ^ = ! , $ = -2tt , 

the analogue of eqs( l4.3ip read like, 

K = , (Co,5 ) = (l,0) , (C M ,5 M ) = (1,0) 

k' = ~ , (cs,5£) = (i,o) , (c;, s;) = (o, 1) 

Putting these values back into (14.291) . we get the following constraint relations, 



(4.38) 



(4.39) 



which can be also put into the form 

4 

£(i$ + iJ£) = (-l + i)flJ. (4.41) 

i/=i 

These condition are manifestly satisfied by the Borigi-Creutz matrix given by eq( ll.2l) 
and the moduli space of the tensor is given by uu^uu) • Similarly, in the case the 
two zero modes are chosen like 

$ 1 = $ 2 = $ 3 = 0, $4 = 0, $o = 

$; = $' 2 = $'3 = 0, $' 4 = 7T, $ = -7T , 

as in the KW fermions, the analogue of eqs( 14.31p read like, 

K = , (Co, So) = (+1,0) , (C 4 ,S 4 ) = (+1,0) 

K' = , (0^) = (-1,0) , (Ci,5i) = (-1,0) ' 
and for fj, — 1,2, 3, we have 

(^,^) = (1,0) , 

(^,^) = (1,0) . 
The constraint relations (I4.29P become, 



(4.42) 



(4.43) 



(4.44) 



^ + E^ = ° > Rl + Rl + R l = K + R i > ( 4 - 45 ) 

exactly as for the ffVF matrix given by eq (11.3j) : see also next section. Notice that, in 
this example, the 4+4 constraint relations are all of them imposed on R l , the real part 
of O.^; so the moduli space associated with these classes of solutions is given by 



Wl> - g ° (4 ' 3) x ^° (4 ' 5) (446) 
32 ~ SO (4) x SO (3) SO (4) x SO (5) 1 } 

Below, we give three claims summarizing the basic results following from this construc- 
tion. 



21 



4.2.2 three claims 



Claim 1 : Moduli space o/X>k,n 

The moduli space of the operator ~ ^ 7 M f^e mk ' Ai + foe) on \D hyperdia- 
mond is given by the complexification of the real 20- dimensional space 9^20 = so(4^^xSO(5) • 
E'oc/i copy o/ the two 9^20 ' s associated with the real R l and imaginary J 1 ^ parts of the 
complex tensor Q l . The combination of these spaces leads to the following complex 20- 
dimensional coset group manifold 

1/(4,5) 



-20 



9^40 (4.47) 



U(A)xU (5) 
Claim 2 : Engineering zero modes 

Each zero mode ofD^Q, located at some wave vector point K = (Ki, K 2 , K 3 , K4) in the 
reciprocal space of the J^D hyperdiamond, requires four real constraint relations given by 
the vanishing conditions Tr [ 7/i 2\n] = 0. These constraints reduce the moduli space 
9Jt4o down to some submanifolds VJI^q. Similarly, two zero modes, located at two points 
K and K' in the reciprocal space, require 4 complex (4 + 4 real) conditions given by 

Tr [7^,0] = , Tr [ 7 ^,n] = . (4.48) 

These constraint relations reduce the moduli space 9Jt 40 down to real 3 2- dimensional 
subspaces 9Jt 32 . In the case of BC fermions, the 9Jl 32 manifold is given by 

U(4,4) 
16 U (4) x U (4) ' 

and in KW fermions; it reads as so S (4)x'so(3) x so^4)x50(5) • ^ or ^ e 9 ener ^ c case °f ( n + m ) 
zero modes ofT>K,ci, one should have 4 (n + m) conditions; and the corresponding moduli 
spaces have the form 9Jt 4 o-4 n _4 m . Typical examples of these submanifolds are given by 

9^20-4n x 9^20-4m With 9^20-4s = SO(^xSO{5-s) ' 

Claim 3 : Continuum limit ofV q+ x,n 

In the neighborhood of each zero mode K with Tr \^^Vk,q\ = 0, the operator T> q+ x,a 
behaves, at first order in the lattice spacing parameter a, like 

4 



(4.49) 



4 

-a Y [R l „ cos aKi - j\ sin aK t ] {qif + O (a 2 ) 



(4.50) 



where the wave vector p^ is related to qi as follows 

4 

V, = E W 5 " + J » S ») 9u + (B%S + J^Co) g , (4.51) 
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with qi = q.A;, K[ = K.A; and C\ = cos aKi, Si = sinaKi. To interpret the q^'s as 
the real wave vector components of the continuum limit, we have to impose moreover 
(R;S U + J»Q) ~ S; and (R^S + J°C ) = 0. 

5 Re-deriving BC and KW fermions 

In this section, we use SU (5) symmetry of the hyperdiamond to re-derive the BC 
fermions and the KW ones from our general action. We first consider BC model and 
then the KW fermions 

5.1 BC fermions 

A simple way to re-derive the BC fermions from our proposal is to work in the reciprocal 
space and compare the Dirac operators following from the models (12.11) and (11.11) . In the 
case of the BC fermions, the Dirac operator T>bc reads, up to a scale factor and upon 
dropping the mass term, as follows 

i 4 2i i 4 
£>bc ~ -^7 M sina^ r V] - V) cos ap^ (5.1) 

with 



a *■ — * a a 



which we rewrite as 



| (7 1 + 7 2 + 7 3 + 7 4 ) , (5-2) 
/ 1 \ 



4 

^1 



^=1 



1 
1 

v 1 / 



(5.3) 



In our proposal (jl.ip . the Dirac operator reads as 



4 4 

-i Y ^ K + n i\ cos a p° - i H ^ H + %] cos a ^ 

fl=l fl,V=l 



(5.4) 



By equating (15. ip and (15. 4p . one gets the tensor Q l that define the BC model. Taking 
the components of the tensor Q l as follows 

^ ^ (5-5) 
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with 



and putting back into ( 15. 411 . we first get 

4 4 
P = 2^ E 7 ^ ~ SillaP ° + 25 E 7 ^ Silla ^ 



fl,U=l 



Then equating with (15. ip . we obtain: 

£ M - In =0 , or sin ap = 



and 



/i 2 M > 



where £ and are still to determine. Solving these quantities as follows, 



where is as in ( 15. 3ft and where 



/ 1 


1 


1 


1 \ 


1 


1 


1 


1 


1 


1 


1 


1 


V i 


1 


1 


1 / 



then putting back into ( 15. 7ft . we end with the following factorized form 
i 



V 



4 



i 7 M sin ap M — i cos ap [ ^ 7 M ^ ) — % ^ 7 M cos ap M 



4 

E 



E 



tr cos ap^ 



V J/=1 



+IE 

Using (15. 3p . we obtain 

4 

£> = ^ ^7 At sinap A1 - ^rcosapo 



i/=l 



cos ap^ — ~ 7^ cos ap M . 
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But this operator is precisely the one given by (I2.20p with Ts = — 2iT and ap^ = apo. 
Comparing with (15. ip . we discover that we should have 

cosap = 1 =>• p = 0, mod (^) . (5-14) 

Let us discuss the meaning of these constraint equations. In our proposal, the compo- 
nents (po,Pfi) form a 5-dimensional vector pi given by 

4 

Pl = J2 k »- X i > Z = 0,...,4 , (5.15) 

where Af are given by eqs (13. 21I3T5]) . Because (Ag + A] 4 + \% + A3 + A4) = 0, these mo- 
menta satisfy the conservation law 

Po = ~ (Pi +P2 +P3 + Pi) • (5.16) 

Now using the fact that the zero modes of the Dirac operator of the BC fermions are 
either p\ = P2 = P3 = Pa = or p\ = P2 = P3 = Pa = -, one finds indeed that po = 0, 
mod^. 



5.2 KW fermions 

In the KW fermions, the Dirac operator T>kw is given by eq (l2.28l) . By dropping the 
bare mass m , it reads, up to a global scale factor, as follows 

4 3 
V kw ~ I ^ ^ sin ak i* ~ « 7 4 S cos ak v + 1*74 • (5-17) 

/1=1 fl=l 

To make contact between this 4x4 matrix operator and the T> one following from our 
proposal, we should solve the equation V KW = V with , 

4 4 

v ~ S^ M K- n °] sin a ^ + 111 ^ K - K] sin a p» 

(5.18) 

By choosing the components as follows, 

K ± % (5-19) 

with 

fi ° = c, , % = % & + \% > ( 5 - 2 °) 
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or equivalently 





/ 2d 




©2 


©3 


©a 


1 


2C 2 


e? 




©1 


©i 


2 


2C 3 


e? 


03 


i + e§ 


©i 




V 2C 4 


©i 


©2* 


* + 0| 


©1/ 



(5.21) 



where C M and 0^ real. Putting back into (I5.18p . we get the reduced operator, 



V ~ ^^T^sina^ - f cosapo $^ 7 %^ 



4 / 4 

Comparing with (I5.17p . it follows that the choice 



/ \ 




\c J 



( 








\ 


























I ox 


c 2 


c 3 


c 4 / 



leads to 



V ~ i$^ 7 M sina ^ -7 4 [f cocosapo] 



"74 



i^c M cosa^ . 



By requiring the KW zeros 



(1) : Pl =p 2 =p 3 = , p 4 =j9 = 

(2) : pi = p 2 = p 3 = , pi = p Q = n 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



obeying ^2^ =0 Pi = mod(27r), we obtain the following constraints on the coefficients 



(1) : Ci + c 2 + c 3 + c 4 + 2c = , 

(2) : ci + c 2 + c 3 - c 4 - 2c = . 

A particular solution of these relations is given by 

ci = c 2 = c , c 3 = c 4 = -2c , c = 1. 



(5.26) 



(5.27) 
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6 Conclusion and comments 



Motivated by studies on lattice QCD and using results on J^D graphene, we have de- 
veloped in this paper a class of J^D lattice fermions that live on the hyperdiamond L 4 
with typical action given by (11.11) . This lattice action involves a complex tensor Q l that 
captures basic properties on hyperdiamond fermions; in particular: 

(1) it links the J^D Euclidean vector /i and the 5-dimensional A; allowing to exhibit ex- 
plicit both the underlying SO (4) and SU (5) symmetries of the lattice action, 

(2) it encodes the data on the zero modes of the Dirac operator T>. The constraint 
relation giving the zero modes of the Dirac operator reads, in terms of the wave vector 
K, as follows: 

Re (y2^e iaKXl \ = (6.1) 
and its solutions depend indeed on Q! and the weight vectors A/, 

(3) it contains as particular examples the BC and KW fermions; the corresponding 
(fijj and (fijj RW were given in the introduction. 

We end this study by noting that one can use the results developed in this paper to 
investigate the renormalization properties of the gauged version of (11.11) along the same 
lines as done in [5j[6]. For instance, the free quark propagator following from our proposal 
is given by 

yF > D 2 -D 2 + (aM f 

with bare mass M ; and 

4 4 

D fl = ^2tt^e iap Xl , = ^ nJ t e -<op - A| , (6.3) 

1=0 1=0 

and where D 2 = D^D^^ D 2 = D^D^. The gauging of (11.11) may be done by making Q 1 ^ 
a local field; i.e by allowing it to depend on the position. Progress in this direction will 
be given in a future occasion. 
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